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?3 . Abstract 

£ . 

Using the generalized Langevin equation formalism and the process of contraction of the descrip- 
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o 
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tion we derive a general memory function equation for the thermal fluctuations of the local density 
of a simple atomic liquid. From the analysis of the long-time limit of this equation, a striking 



equivalence is suggested between the long-time dynamics of the atomic liquid and the dynamics 
> ■ 

of the corresponding Brownian liquid. This dynamic equivalence is confirmed here by comparing 



molecular and Brownian dynamics simulations of the self-intermediate scattering function and the 



long-time self-diffusion coefficient for the hard-sphere liquid. 
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I. INTRODUCTION. 

In many respects, the properties of colloidal fluids resemble almost perfectly those of the 
correspondent atomic liquid l|-|4| ■ It is well known that the equilibrium phase diagram, and 
in general all the equilibrium thermodynamic properties, of a specific model system (say a 
Lennard- Jones liquid) will be independent of the microscopic (either molecular or Brownian) 
dynamics that govern the motion of the N interacting particles that constitute the system. 
This implies that these equilibrium properties can be generated using either molecular or 
Brownian dynamics simulations [5j. Furthermore, although time-dependent properties are 
expected in general to depend on the specific microscopic dynamics, some features associated 
with the long-time dynamic behavior of the system also seem to be rather insensitive to the 
microscopic short-time dynamics. This appears to be particularly true regarding the rather 
complex dynamic behavior of these systems as they approach the glass transition p-l8]. 
Determining the range of validity of this analogy continues to be a relevant topic in the 
study of the dynamics of liquids. 

From the theoretical side one would like to unify colloidal and atomic liquids in a common 
theoretical description of the relaxation dynamics of the local density fluctuations, which ex- 
plicitly exhibits the origin of the similarities and differences in their macroscopic dynamics. 
One possible general framework for such theoretical analysis is the concept of the generalized 
Langevin equation (GLE) [9|, [lOJ. This equation describes the dynamics of the thermal fluc- 
tuations 8a,i(t) (= Oj(t) — a^ q ) of the instantaneous value of the macroscopic variables a,i(t) 
(i = 1, 2, ..., z/), around its equilibrium value a^ 9 , and has the structure of the most general 
linear stochastic equation with additive noise for the vector Sa(t) = [5di(t), 5a,2(t), ..., 5a u (t)] 
(with the dagger meaning transpose). The GLE equation has been widely used in the de- 
scription of thermal fluctuation phenomena in simple liquid systems, and Boon and Yip's 



textbook JjJ contains a detailed account of its early use to describe the dynamics of sim- 



ple liquids. Although this stochastic equation is conventionally associated with the Mori- 
Zwanzig projection operator formalism [12j, ll3(, in reality its structure is not a consequence 
of the hamiltonian basis of Mori-Zwanzig's derivation; instead, it is essentially equivalent to 
the mathematical condition of stationarity [9|. 

Thus, in Ref. [14| the GLE formalism, understood in the latter manner, was employed 
to derive the most general diffusion equation of a model Brownian liquid (i.e., an idealized 



monodisperse colloidal suspension in the absence of hydrodynamic interactions) formed by 
N spherical Brownian particles interacting between them through direct (i.e., conservative) 
forces, but in the absence of hydrodynamic interactions. The resulting general memory 
function expression for the intermediate scattering function (ISF) F(k,t) and for its self 



component Fs(k,t), were later employed in the construction o 



alized Langevin equation (SCGLE) theory of colloic 



dynamics 



the self-consistent gener- 
ly applied 
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21] to the 



to the description of dynamic arrest phenomena 17H19| and more recently 20|, 
construction of a first-principles theory of equilibration and aging of colloidal glass-forming 
liquids. 

With the aim of investigating the relationship between the dynamics of atomic and Brow- 
nian liquids, here we start the extension of these theoretical developments to describe the 
macroscopic dynamics of both kinds of systems within the same theoretical formalism. With 
this general intention in mind, in the present paper we discuss the application of the gen- 
eralized Langevin equation formalism above, to the derivation of general memory-function 
expressions for the (collective and self) intermediate scattering functions of an atomic liq- 
uid. These expressions should in principle be capable of describing the crossover behavior 
of these properties between their ballistic short time limit and their diffusive long-time be- 
havior. Although in practice we do not use these expressions here to numerically evaluate 
these functions in the short- or intermediate time-regime t ~ r (where To is the mean free 
time), we find that in their long-time limit, t ^> tq, these expressions for F(k, t) and F$(k, t) 
become essentially identical to the corresponding expressions for a colloidal fluid, strongly 
suggesting a well defined long-time dynamic correspondence between atomic and colloidal 
liquids. 

The strategy that we shall employ to derive the memory function equations for the 
intermediate scattering functions of our model atomic liquid will actually rely very heavily on 
the referred previous derivation [14J of the time-evolution equations for F(k,t) and Fs(k,t) 
of the corresponding idealized Brownian fluid. The rationale for this is the rather simple 
observation that the essential difference between an atomic liquid and its idealized Brownian 
counterpart (a colloidal liquid in the absence of hydrodynamic interactions) is the presence, 
in the microscopic equations of motion of the latter, of the friction force — C Vj(t) due to 
the supporting solvent and the corresponding fluctuating force fw(£). Thus, we first review 
the derivation of Ref. [14J . with the aim of keeping track of the effects of these friction 



terms. This aspect of the present work is developed in section [TTJ At the end of the section, 
we simply take the C — > limit of the end result of the referred derivation, to obtain the 
corresponding time-evolution equations for F(k, t) and F s (k, t) of our atomic liquid (namely, 
Eqs. flPgD and (JX30]) ). 

The next task of this work is to analyze the long-time limit of these results for F(k, t) 
and Fs(k,t). In Ref. 14|, dealing with Brownian systems, this limit was referred to as the 
"overdamped" limit, corresponding to times t much longer than the relaxation time t^ = 
M/(^ of the velocity autocorrelation function. This relaxation results from the damping of 
the particle's momentum due to the friction force — C Vi(£). Thus, in that case t^ sets the 
crossover timescale from the early initial regime t « t^ s \ where the inertial effects are still 
important, to the long-time regime t » r^ s \ where the motion of the suspended particles 
is purely diffusive, and described by the short-time self-diffusion coefficient D^ = ksT '/C • 
In contrast, in atomic liquids an analogous timescale is apparently absent, since there is 
not any material solvent exerting damping friction forces. In spite of that, in Section IHH 
we analyze the long-time limit of the time-evolution equations for F(k,t) and Fs(k,t) of 
the atomic liquid derived in Section [Til We find that in this limit, these equations happen 
to adopt the same structure as the corresponding equations for Brownian systems in their 
overdamped limit. As a result of this analysis, we conclude that the parameter playing the 
role of the short-time self-diffusion coefficient D^ is now the self-diffusion coefficient D° 
determined by kinetic theory. 

This formal dynamic correspondence has important physical consequences, expressed in 
terms of well defined scaling properties of the dynamics of two fluid systems which only differ 
in the microscopic laws that govern the motion of the constituent particles (either molecular 
or Brownian dynamics) . The most relevant of such consequences are briefly discussed in the 
final section (Section [V} of this paper. 

II. ATOMIC FLUID AS A FRICTIONLESS BROWNIAN LIQUID. 

Let us start by reviewing the derivation in Ref. [14j of the time-evolution equations of 
F(k,t) and Fs(k,t) of an idealized monodisperse colloidal suspension in the absence of hy- 
drodynamic interactions, formed by N spherical particles in a volume V, whose microscopic 



dynamics is described by the iV-particle Langevin equations 22M24 ] 



M dv^)_ ^ _ cis)vi{t) + f (s ){t) + j2 Fij{tl {i = 1,2,. ..,N). (2.1) 

In these equations, M is the mass and Vj(t) the velocity of the zth particle, and Q 8 ' is its 
friction coefficient in the absence of interactions. Also, fj (£) is a random force, modeled as a 
Gaussian white noise of zero mean, and variance given by (fv(£)f ■(())) = ksT (,^25^)8^ I 
(z, j = 1, 2, . . . , JV; I being the 3x3 unit tensor). The direct interactions between the par- 
ticles are represented by the sum of the pairwise forces Fy that the jth particle exerts on 
particle i, i.e., Fy is obtained from the pair potential w(|i*j — r^- 1 ) . 

Our goal is to derive the macroscopic time-evolution equations for the ISFs F(k, t) and 
F$(k,t), starting from this microscopic level of description. Some of the most important 
features of such general time evolution equations for F(k,t) and Fg(k,t) can be written, 
however, right at the outset, since they derive from the general selection rules [9J originating 
from the stationarity condition and from other symmetry properties of the macroscopic vari- 
ables whose dynamics couple to the dynamics of the local particle concentration. This was 
the approach adopted in Ref. [l4j, which derived the most general time-evolution equation 
for the fluctuations of the local concentration n(r, t) of colloidal particles, consistent with 
the selection rules referred to above. The specific information of the microscopic dynamics, 
was then employed in the approximate or partial determination of those elements of the 
time-evolution equation that such selection rules left undetermined. This section briefly 
summarizes the main steps of such derivation. 

At each step of the following derivation, however, we urge the reader to keep track of 
the particular case in which the friction term — £( s 'Vi(t) and its corresponding fluctuating 
force fj (£) are absent, and to recognize that an atomic liquid can be viewed as the present 
Brownian liquid in the limit of an infinitely tenuous solvent, such that the Stokes friction 
coefficient (^ vanishes. Thus, we shall take the limit £^ — > in the general memory- 
function expressions for F(k,t) and F$(k,t) derived in this section. In such limit one is left 
only with the particles in the vacuum, and these equations for F(k, t) and F$(k, t) will then 
become the exact memory function expressions for the ISFs of an atomic liquid. 

Thus, let us first recall that the basis of the GLE formalism are the general mathematical 
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conditions stated by the theorem of stationarity |9|. This theorem states that the equation 
describing the dynamics of the thermal fluctuations Sai(t) (= di(t) — a^ q ) of the instantaneous 



value of the macroscopic variables cti(t) (i = 1, 2, ..., v) around its equilibrium value a^ q must 
have the structure of the most general linear stochastic equation with additive noise for the 
vector 8&(t) = [5a,i(t), 5a 2 (t), ..., 5a u {t)} , namely, 



dSa(t) 



X 

dt -u X - l 8*{t) - J L(t - t'^Sait'W + f (£). (2.2) 

o 

In this equation x is the matrix of static correlations, Xij — (^ a i(0)5a*(0)), ui is an anti- 
Hermitian matrix (uij = —00%), and the matrix L(t) is determined by the fluctuation- 
dissipation relation L^t) = (fi(t)fj{0)), where fi(t) is the ith component of the vector of 
random forces f(t). 

For the present purpose, we choose the components of the state vector 8&(t) as 

8a(t) = [5n(k,t),<5j(k,t),(5^(k,t),^(k,t)] t , (2.3) 

with the following definitions. First, di(t) is the Fourier transform Sn(h,t) of the fluctua- 
tions Sn(r, t) = n(r,t) — n of the local concentration n(r,t) around its bulk value n. The 
microscopic definition of <5n(k, t) (for k = 0) is 

1 N 
Sn(k,t) = ^=E eik ' riW ' ( 2 - 4 ) 



N 1 

where rj(t) is the position of the ith colloidal particle at time t. Normalized in this manner 
5n(k,t) is such that its static correlation is Xnn(k) = (Sn(k, 0)Sn(— k, 0)) = S(k), where 
S(k) is the static structure factor of the bulk suspension. 

Taking the time-derivative of <5n(k, t) we have the continuity equation, 

dSn(k, t) 



dt 



ikSjtfat), (2.5) 



where Sji(k,t) = ji(k,t) = k-j(k, t) is the component of the current j(k, t) in the direction 

k of the vector k, i.e., 

1 N ~ 
3l (k,t) = ^=J2^Mt)e^ (t) (2.6) 



N . . 



with Vi(t) = d,Ti(t)/dt. Thus, d2(t) = Sji(k,t), whose static correlation matrix is 

X 33 = k B T/M. (2.7) 



If we take the time-derivative of the current in Eq. f |2.6|) . and employ the iV-particle 
Langevin equation, Eq. (12.11) . we are led to the following result 

N -^ ,. s ., N 



2 

g*T»(t) 



at m -" v ' ; m v 

v ~ i=i ■ ~ i=i 

(2.8) 

where 

1 * 
/ w (k,*) = ^E £f 5 (*) e<k " r ' (t) - ( 2 - 9 ) 

and Fj(t) = X^jFij(t). This equation can also be written as 

^ = -— 5j z (k,t) + — — — +jfcfo (k,£), (2.10) 

with 5a zz (k,t) being the instantaneous fluctuation of the isotropic diagonal component of 

the stress tensor 

1 N ( 1 a ? ^\ 

""(M - ^ E W - ^ E ^^)| e^ W , (2.H) 

where 

In these equations, r^ = r^ — r,,-, and ufaj) is the pair potential. 
Let us now write 5a zz (k, t) as 

8a zz (k,t) = 8p(k,t) + Sa K (k,t) + 8au(k,t), (2.13) 

with Sp(k,t) = [xjj/S(k)]5n(k,t) being the Fourier transform of the local pressure fluctua- 
tions, and with Sctk (k,t) and Sau(k,t) being the statically orthogonal kinetic and configu- 
rational components of [5cr zz (k,t) — Sp(k, t)], defined as 

1 - 
Sa K (k,t) e-J^V^^O-^XM), (2.14) 



N . 

4=1 

and 

1 N r a r P 

5au(k,t) = -— _ J2Y, -^Pkir^e***® - Sp(k,t) + XjJ Sn(k,t). (2.15) 

* j=l j^i l J 

This completes the microscopic definition of the components Sn(k, t), Sj(k, t), 8(7k (k, t), and 
Scru{k, t) of the state vector a(£), which are then found in Eqs. (12.4[) . (12. 6p . (12.14[) . and 
(I2.15p . respectively. 



As a result, we finally rewrite the momentum conservation equation, Eq. (I2.10p . as 
d5j t (k,t) 



Of 



^rta(M) + ]^/ (s) (k,t) +ikSp(k,t) +ik5a K (k,t) +tk5au(k,t). (2.16) 



This equation, together with the continuity equation in Eq. ( 12. 5p . couple the variables 
5n(k,t) and Sj(k, t) with the variables <5crx(k, t) and 5ajj(k., t), whose time-evolution equa- 
tion must now be determined, and the GLE formalism provides a natural manner to do 
that. For this, one first performs a straightforward statistical thermodynamical calculation 
of the matrix x °f static correlations Xij — ( a j(0)c^(0)), with the following result [14J 



X 



with Xr. 



Un 

xa ° ° 

xkk 
xvu 

S(k) and Xjj — ksT/M, and with xkk and xuu given by 



(2.17) 



Xkk 



2X?, 



and 



Xuu 



x 2 - 



1 + n I drg(r) 



d 2 /3u(r) fl — cos(kz) 



(2.18) 



(2.19) 



dz 2 \ k 2 J S(k)_ 
We then write up the generalized Langevin equation for our vector Sa(t) in the format of 

Eq. ( 12. 2p . For this, we first notice that all the variables, except 5a,2{t) = 5ji(k,t), are even 

functions under time-reversal. According to Onsager reciprocity relations, and the general 

anti-hermiticity of to and hermiticity of L(z) [9], we have that the only possibly non-zero 

elements of the matrix u and Liz) are 



to- 



L(t) 






U nj 





-<J 


UJjK 


LOju 


- 


-u* K 





- 


-sv ° 





J- J nn 


L n K 


L n U 





Ljj 





T* 

^nK 


L KK 


Lru 


T* 

. ^nU 


L* KU 


Luu _ 



(2.20) 



(2.21) 



The determination of the non-zero elements of u and of some of the non-zero elements of 
L(t) is rather straightforward, since, from the exact continuity equation, 

^f^ = ik Sjl (Kt) (2.22) 

we immediately see that u n j = —ikxjj, and that L nn = L nK = L nU = 0. Similarly, from 
eq. ( 12.161) we can see that oJjkxJu = u juXuu = ~^ anc ^ ^jjXjj = C^V^- As a result, 
all the elements of the "frecuency" matrix u have been determined, and in fact, only the 
kinetic coefficients L KK (k,z), L KU (k,z) = L UK (k,z), and L V u{k,z) remain undetermined 
by general symmetry principles, or physical principles such as mass or momentum conserva- 
tion. Thus, the time-evolution equations that complete the non-contracted description for 
the components of the vector 5a(t) are the mass and momentum conservation equations, 
Eqs. ( I2.5P and ( 12.161) . along with the time-evolution equations for Saxik, t) and 5au{k, t), 
namely, 

cMcr^k, t) 



dt 



ikxKKXjj^JiiK t) - / L KK (k, t - t')x K 1 K S(r K (k, t)dt' 
Jo 

- [ L UK (k, t - t')xuu^u{^ t)dt' + f K (k, t) (2.23) 

Jo 



and 



dSau (k, t) 
dt 



ikxuuXjj~$3i(k,t) - / L uu (k,t-t')x u 1 u Sa u (k,t)dt' 
Jo 

- [ L UK (k,t- t')xK K 5a K (k,t)dt' + fu(k,t). (2.24) 

Jo 

In these equations, only LxKik-, t), Luu(k,t), and Lux{k, t) remain unknown. 

The extended dynamic description provided by Eqs. (12.51) . (I2.16p . (I2.23p . and (I2.24p can 

now be contracted down to a single time-evolution equation involving only Sn(k, t) (9J. This 

essentially amounts to formally eliminating the variables 5j(k,t), 5<7K{k,t), and 5au{k, t), 

from this system of equations. The result of such contraction procedure reads [14| 



05n(k, t) 



f L(k,t -t')XnnHKt')dt' + f(k,t), (2.25) 

Jo 



dt 

where f(k, t) is a random term with zero mean and time-dependent correlation function 
(/(k, 0)/(— k, 0)) = L(k,t) with L(k,t) given, in Laplace space, by 

L(k, z) = r— ^ ; -. (2.26) 
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with z^ = C^/M and 

ALjj(k, z) 



k 2 XKK 



k 2 Xuu 



+ 



-KUXyu 



z+L 



+LkkX K k 



z + LkkxA z + Luuxuh 



X K K L KU L UKX UU 



(2.27) 



z+LkkXkk 

Multiplying Eq. (I2.25P by 8n{— k, 0), and taking the equilibrium average, this equation 
becomes the time-evolution equation for the intermediate scattering function F(k, t) = 
(Sn(k, t)Sn(— k, 0)), an equation that can be written as an expression for the Laplace trans- 
form F(k,z) in terms of the memory functions L KK (k,z), L KU (k,z) = L UK (k,z), and 
L uv (k,z), namely, 

S(k) 



F(k,z) 



z + 



k^s-Hk)xjj 



(2.28) 



Z+Z (s )+ **» XKK 



k 2 X j:j XUU 



2 + L KKX K K 



-KU^uu 



Z + L KK* KK 



i r -1 X KK L KU L UK*Tjn 

Z + L KKX KK 



At this point we can discuss the limit of vanishing solvent friction, ^ — > 0. As discussed 
above, in this limit our Brownian fluid becomes a Newtonian system, in the sense that its 
microscopic dynamics is described by Eq. (12.11) without the friction and fluctuating terms. 
Thus, the expression for F(k, z) describing the collective dynamics of an atomic liquid can 
be obtained from the previous expression by simply setting z^ = 0, i.e, 

S(k) 



F(k,z) 



tfS-HVxn 



(2.29) 



z+- 



k z 



X jj X KK 
z + l KKXkK 



k2 x^xuu 



^KUXjju 



Z + L KKX KK 



z + L uu*uu ~r^ — -1 

z+ L KKX KK 



In a completely analogous manner we can derive the corresponding expression for the 
self-lSF Fg(k,t), with the following result 



F s (k,z) 



1 



z + 



'Xjj 



(2.30) 



^ 2 A 



-1 (S) 
3 X UU 



•(S) (S)-l 
'Kl' x UU 



Z + L KK*l 



(S) 



with 



X.uu 



z+L 



n x% 

k 2 



(S) (S)-l_ 



1 (S) (S) (S)-l 
'Jtg KU UK X UU 



*-+l¥Ix 



1 



KK^KK 



drg(r) 



d 2 (5u{r) 
dz 2 



(2.31) 



These general results now will serve as the basis for the analysis of the long-time dynamics 
of an atomic liquid, carried out in the following section. 
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III. LONG-TIME DYNAMIC EQUIVALENCE OF ATOMIC AND COLLOIDAL 
LIQUIDS. 

In this section we analyze the long-time (or small frequency) limit of the general ex- 
pressions for F(k,z) and F s (k,z) in Eqs. ( I2.29P and (I2.30p . With this purpose, as an 



additional approximation (following Ref. [14|, but introduced here only for simplicity) let 
us first neglect the possible crossed kinetic couplings represented by the memory functions 
Lxu{k, z) = Lux{k, z) in this equation. This leads to simpler expression for the ISF of an 
atomic liquid, namely, 

W*, z) = #L- (3-1) 



z + 



— 9 — ^"T ~? — ^"T 

k X:a XKK , k Xji *UU 



=1 — i =r 

z + L KK( k ' z )x KK z +L UU (k,z)x uu 



and 



F s (k, z) = s . (3.2) 

~ k2 X 3 ^XKK | fc%/x^ 

Eqs. ( 13. ip and ( 13.21) express F(k,t) and Fs(k,t) in terms of the unknown mem- 
ory functions LxK(k,z), Luu{k,z), L KK (k,z) and L v {j{k,z). To understand the prop- 
erties of these memory functions, with the aim of introducing additional approximations 
or simplifications, it helps to analyze their physical meaning. For this, let us recall 
that the memory functions Lxx(k,z) and L KK (k,z) are associated with the relaxation 
of the kinetic part cr^f(k, t) = iV -1 / 2 Yli=i v"v^ e t]i ' ri ^ of the stress tensor, whose trace 
crR-(k, t) = N~ 1 / 2 Y2i=i v 2 e jk ' r ^*- ) is directly related with the FT of the local kinetic energy 
density. Thus, LKx(k,z) and L KK (k, z) essentially describe the transport of molecular ki- 
netic energy, i.e., the transport of heat. These transport processes occur primarily by means 
of molecular collisions and quickly lead to a uniform distribution of the mean kinetic en- 
ergy of the particles, i.e., to thermal (but not thermodynamic!) equilibrium. As a result, 
these memory functions may be expected to be related with heat conductivity, and to decay 
within molecular collision times. The memory functions Luu(k,z) and L\jjj(k,z), on the 
other hand, describe the relaxation of the configurational component of the stress tensor, 
which involves structural relaxation processes that may decay after much longer relaxation 
times. 

Because of this, if one is interested in the long-time behavior of the ISFs, one may neglect 
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the frequency- dependence of LrKK(k,z), and replace it by its zero-frequency limit, 

L KK (k, z) « L K (k) = lim L KK (k, z) 
in Eq. (13. ip . and similarly for L KK (k,z), 

Lf K {k,z) « 4 5) (A;) ee limL<&(fc,z), 



(3.3) 



(3.4) 



in Eq. (I3.2p . In addition, we also assume that the kinetic coefficients L^(k) and L K (k) are 
not fundamentally different from each other, so that we neglect their possible differences, 



L K (k)*iL$>(k). 



(s). 



(3.5) 



At this point we take the desired long-time limit t » To in the resulting approximate 
expressions for F(k,z) and Fs(k,z). This amounts to neglecting the frequency z compared 
with the frequencies zd = LK K (k,z)x 7 K 1 K and z B = ^xJ^Xkk/zd in Eqs. (13.11) and (13.21) . 
which leads to the "overdamped" form of these expressions, namely, 

S(k) 



and 



F(k,z) 



F s (k,z) 



, fc 2 5~ 1 (fc)J' 
Z + 1+C(k,z) 



Z + 



k 2 D° 



(3.6) 



(3.7) 



1+C s (k,z) 

where we have defined the memory functions C(k, z) and Cs(k, z) as 



and 



C(k,z) = 



C s (k,z) = 



k 2 D X T 2 xuu 



z + Luu(k,z)xuu 



k 2 D\-fx^u 



z + ^uu\k) z )Xuu 
respectively. 

In these equations we have denoted the unknown frequency z D = L^ K '{k)xxK as 



(3.8) 



(3.9) 



L { * ) (k) X - K 1 K = 2k 2 D°. 



(3.10) 



The use of the symbol D° is, of course, not accidental, since this parameter can be identified 
with the self- diffusion coefficient that describes the sequence of ballistic random flights of 
a tracer particle as it collides with its neighbor particles. To see this, notice that in the 
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conditions in which the effects of the configurational memory function Cs{k, z) are negligible 
(such as in the low-density regime, in which Xuu = Xuu = 0), Eq. (13.71) becomes 

F »<*-*> M 7+k5' (311) 

or 

F s (k,t)^e- k2D0t . (3.12) 

This result implies that the MSD is given by W{t) ~ D°t, i.e., that the motion of a tracer 
particle after many collision times will be diffusive. The corresponding diffusion coefficient 
D° must then be identical to that determined by kinetic-theoretical arguments, i.e., must be 
given by D° = (l ) 2 /r , where Iq and Tq are, respectively, the mean free path and the mean 
free time. Since Zo/ r o = v 0i D° can also be written as D° = v l . If we then estimate the 
mean free path Iq to be given by Iq ~ 1/na 2 , with n = N/V and with a being the collision 



diameter of the particles, we then have that D° ~ y^ksT/M /(na 2 ). In fact, the rigorous 



value of D° is 



25] 



n 3 fk B T\ 1/2 1 

D = —^ -V — o- (3.13) 

80F V M ) na 2 K J 

The comparison of the overdamped expressions for F(k, z) and F$(k, z) in Eqs. ( 13. 6^ - 
( 13. 9p above, with the corresponding overdamped results of a colloidal liquid (i.e., with Eqs. 
(4.24) and (4.33) of Ref. [14J), reveals the remarkable formal identity between the long- 
time expressions for F(k,t) and F$(k,t) of an atomic liquid, and the corresponding results 
for the analogous colloidal system. The fundamental difference between these two cases is 
to be found in the definition of the diffusion coefficient D°, which in the present (atomic) 
case depends on temperature and density, and is given by the kinetic-theoretical result in 
Eq. ( I3.13p . whereas in colloidal liquids it is a constant, identical to the short-time self- 
diffusion coefficient given, for example, by the Einstein-Stokes expression in the absence of 
hydrodynamic interactions. Thus, this formal identity implies that the long-time dynamic 
properties of an atomic liquid will then coincide with the corresponding properties of a 
colloidal system with the same S(k), provided that the time is scaled as D°t, with the 
respective meaning and definition of D°. This observation has important implications, which 
can be tested, for example, by comparing the simulation results for Fs(k, t) obtained by both, 
molecular dynamics and Brownian dynamics, for the same system and conditions. 
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IV. TEST OF THE PREDICTED LONG-TIME DYNAMIC EQUIVALENCE. 

In this section we perform the test of the predicted long-time dynamic equivalence between 
a model atomic liquid and its corresponding Brownian fluid. This dynamic equivalence is 
tested here by comparing the macroscopic dynamics of the hard sphere liquid when the 
motion of its constituent particles is described, respectively, by Eqs. (12. ip without and with 
the solvent friction terms present, i.e., by performing and comparing the molecular and the 
Brownian dynamics simulations of these properties. 

As a reference let us first recall the exact short-time limit of the self-ISF of an atomic 
liquid. Since for correlation times t shorter than the mean free time Tq all the particles move 
ballistically, [r«(t)-r 4 (0)] = v% we have that F s (k, t) « (l/iV)(Eii ex P [« k ' v ;( )*])- Usin S 
the equilibrium distribution of the initial velocities Vj(0), one can see that the exact short- 
time limit of the self-ISF is given by Fs(k,t) = exp(— |fc 2 f 2 £ 2 ). This expression provides an 
excellent approximation at small volume fractions, where Fs(k,t) has decayed to negligible 
values for t ~ r°, as illustrated by its comparison in the main panel of Fig. (CQ) with the 
MD-simulated Fs(k,t) for the hard sphere fluid at <fi = 0.1 and ka = 7.1. 

The MD simulations were conducted on a soft-sphere system, and the results were then 



mapped onto those of the equivalent hard-sphere liquid as discussed in Ref. [26] . The 
soft-sphere simulations were carried out using the velo city- ver let algorithm with N = 1000 
particles of the same mass M in a volume V and a time step At/t* = lXlO~ 3 ^ma 2 /e. 
During the equilibration and production cycles, temperature was kept constant by a simple 
rescaling of the velocities of the particles every 100 time steps. For high volume fractions we 
used polydisperse systems, where the diameters of the N particles were evenly distributed 
between a(l — w/2) and cf(l + w/2), with a being the mean diameter. We consider the case 



w = 0.3, corresponding to a polydispersity s a = w/y/12 = 0.0866. The length, mass, and 



time units employed are, respectively, a, M, and a^M/ksT. The simulations are carried 
out for an array of volume fractions <fi = (ir/6)na 3 where <r 3 is the third moment of the size 
distribution and n is the total number density n = N/V . 

Defining the relaxation time r a by the condition Fs(k,T a ) = 1/e, we have that in the 
ballistic regime r a can be approximated by r a = j^-, which is the low-density limiting value 
represented in the inset of Fig. [1] by the horizontal dashed line. The inset also plots the 
simulation results for r Q in a wide range of volume fractions, to show the deviations from 
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this limiting behavior as the density is increased. Beyond this low-density regime, these 
deviations become increasingly more important, as also illustrated in the main panel of Fig. 
Q]by the MD simulation results for F$(k, t) at the near-freezing volume fractions <fi = 0.4 and 
0.5. Here, of course, the ^-independent limit Fg(k,t) = exp(—^k 2 v 2 t 2 ) is clearly inadequate, 
although the Gaussian approximation, Fs(k,t) w exp[— k 2 W(t)} still provides an accurate 
representation of the short-time decay of this function. This is illustrated by the solid lines 
of the main panel of Fig. (TJ which result from employing the MD-simulated data for the 
mean squared displacement W(t) in Fs(k,t) = exp[— k 2 W(t)}. 




FIG. 1: Molecular dynamics results for the self-ISF Fs(k,t) of a hard-sphere fluid as a function of 
time t (expressed in "molecular" units [c/^o]) a t fixed wave- vector ka = 7.1, and volume fractions 
(j) = 0.1 (solid circles), 0.4 (empty circles), and 0.5 (striped circles). The dashed line is the exact 
limit Fs(k,t) = exp(— ^^v^t 2 ) and the solid lines are the results of the Gaussian approximation 
Fs(k,t) ~ exp[— k 2 W(t)], with W(t) given by the same molecular dynamics simulation data. 
In the inset we plot the relaxation time r a (also in units of [c/uo]), defined by the condition 
Fs(k,T a ) = 1/e, for these and other volume fractions; the horizontal dashed line indicates the 
limiting value VQT a /a = t§ = 0.199. 



With this low-density short-time ballistic limiting behavior as a reference, let us now com- 
pare the simulation results for F$(k, t) obtained by both, molecular dynamics and Brownian 
dynamics simulations, for the same system and conditions. For this comparison, in addition 
to the molecular dynamics simulations, we performed Brownian dynamics simulations of 
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the hard sphere 



QE 



iquid using the conventional Ermak and McCammon's Brownian dynamics 



algorithm |5|, |27J on a soft sphere fluid, and then mapping the results onto those of the 
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. The 



hard-sphere liquid according to the methodology proposed and explained in Ref. 
resulting comparison provides a test of the theoretical prediction of the previous sections, 
that the dynamics of an atomic liquid coincides with the dynamics of the corresponding 
Brownian fluid in the opposite regime, i.e., for high densities and long times. Thus, Fig. 
Et^a), presents both simulation results for the hard sphere system at three volume fractions, 
= 0.50, 0.548, and 0.571, representing the metastable regime of the hard sphere liquid. As 
this figure illustrates, plotting Fs(k,t) as a function of the scaled time t* = D°t/a 2 clearly 
exhibits the expected long-time dynamic equivalence between atomic and Brownian liquids. 
We notice, however, that this long-time dynamic equivalence is not observed in Fg(k,t) at 
lower volume fractions, corresponding to the stable fluid regime (<f> ;$ 0.45). The reason for 
this is that in such regime, illustrated in Fig. [Q the decay of F s (k, t) to a value ps e _1 occurs 
within times comparable to the mean free time tq and is, hence, intrinsically ballistic. It is 
only at higher volume fractions that this long-time dynamic equivalence is fully exhibited 
by the diffusive decay of Fs(k,t), as illustrated by Fig. E(a). 

Another manner to summarize this observation is to compare the volume fraction de- 
pendence of the relaxation time r a of both, molecular and Brownian dynamics. In Fig. 
Hfb) these simulation results are presented in terms of the dimensionless a-relaxation time 
t* = k 2 D°T a . For a Brownian liquid r a — > l/k 2 D° as <f> — > 0, with a ^-independent short- 
time diffusion coefficient D°, so that r* — > 1 as (f> — > 0. As discussed in the previous 
section, however, for atomic liquids r a — > \2/kvo as (f> — > 0, so that in the same limit 
r* — > (ko~)y2ir /16(f), where we have taken into account the fact that in this case, the short- 
time diffusion coefficient D° is given by the kinetic-theoretical result in Eq. ( 13. X3j) . This 
limiting behavior was represented by the horizontal dashed line of Fig. (TJ and is now rep- 
resented by the dashed curve of Fig. |2]^b). From the comparison in this figure one can see 
that the long-time dynamic equivalence manifests itself in the collapse of the molecular and 
Brownian dynamics data for r* at high volume fractions. For smaller volume fractions, the 
differences in the short-time behavior of F$(k, t) lead to the observed differences between the 
molecular and Brownian dynamics results for r* below a crossover volume fraction located 
near the freezing transition of the HS liquid. 

The solid curve in Fig. MJd) is the prediction for r* = k 2 D°r a of the self-consistent gen- 
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FIG. 2: (a) Molecular dynamics (solid symbols) and Brownian dynamics (empty symbols) simula- 
tion results for the self- intermediate scattering function Fg(k, t) of the hard sphere liquid at volume 
fraction (j) = 0.50, 0.548, and 0.571, evaluated at the main peak of the static structure factor and 
plotted as a function of the dimensionless time t* = D°t/a 2 . (b) Volume fraction dependence 
of the dimensionless a-relaxation time r* = k 2 D°T a of the hard sphere liquid determined from 
the corresponding molecular dynamics (solid symbols) and Brownian dynamics (empty symbols) 
simulations. The dashed curve represents the low-density limit r* = (/cc)v / 27r/160, whereas the 
solid curve correspond to the results of the SCGLE theory (Eqs. (l),(2),(5)-(8) of Ref. [19(, with 
k c = 1.305(2tt/<t)). 



eralized Langevin equation (SCGLE) theory of colloid dynamics, i.e., of Eqs. (l),(2),(5)-(8) 
of Ref. |19[. These are actually Eqs. (13. 6p and (13. 7p above, complemented by the closure 
relation C(k,t) = Cs{k,t) = X(k)A((t), where A((t) is the time-dependent friction function 
describing the configurational contribution to the friction force on a tracer particle (given by 
Eq. (6) of Ref. [19]). The static structure factor of the hard sphere system, needed as an in- 
put in these equations, is provided by the Percus-Yevick approximation with its Verlet-Weis 
correction [29], l30t] . The function X(k) = 1/[1 + (k/k c ) 2 ] is a phenomenological "interpolat- 
ing" function, with the cutoff wave-vector k c used here to calibrate the SCGLE theory by 
optimizing the overall agreement of its predictions with the data for r* constituted by the 
totality of the Brownian dynamics results (squares) and by the molecular dynamics data 
corresponding to the metastable liquid (0.5 ^$ <p) in this figure. This calibration procedure 
results in the value k c = 1.305(27r/er). 
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As said above, the short-time differences between the molecular and the Brownian 
dynamics data for r* in Fig. [2(b) appear at densities below a crossover volume frac- 
tion located, for the data in this figure, near the freezing transition of the HS liquid. 
The location of this crossover depends, however, on the wave-vector k at which the de- 
cay of Fg(k,t) is being observed, moving to a vanishing value in the long- wavelength 
limit, k — > 0. This means that in this limit the molecular and Brownian dynamics re- 
sults for r* will be identical at all volume fractions. In fact, this is also what hap- 
pens to the most representative long-time dynamic property, namely, the long-time self- 
diffusion coefficient. D L is defined as D L = lim t ^ 00 ((Ar(t)) 2 )/6t, but is also given by 
D L = lim k ^ \im z ^ [k 2 F s (k, z)}- 1 = D°/[l + C s (k = 0, z = 0)]. According to Eq. (I577J) 
above, and within the SCGLE closure Cs{k,t) = X(k)A((t), for an atomic system this 
parameter, scaled as D* = Dl/D°, can be written as 

poo 

D* = 1/[1 + / AC(t)dt], (4.1) 

Jo 



with A(*(t) given, according to Eq. (6) of Ref. [19|, by 

^w = ^-/* rw> " 11 

3 (27r) n J 



F(k,t)F s (k,t). (4.2) 



S(k) 

These equations, however, are identical to their colloidal counterpart. Thus, they imply 
that the parameter D* of an atomic liquid must be indistinguishable from the corresponding 
parameter of the equivalent colloidal system with the same interactions and the same static 
structure factor. 

The accuracy of this important and distinct prediction can also be checked by comparing 
the corresponding molecular and Brownian dynamics results. Thus, in Fig. [3] we plot 
molecular dynamics data for Dl{4>) of a hard-sphere fluid both, in the "usual" atomic units 
a{k B TjMfl 2 , and scaled as D*(<p) = D L ((j))/D (<p), with D (<j>) given by Eq. (13131) . The 
same figure also presents available Brownian dynamics simulation results for Dl(4>) of the 
hard sphere system without hydrodynamic interactions, also scaled as D*((j>) = Dl(4>)/D°, 
but with D° being the ^-independent short-time self- diffusion coefficient of the Brownian 
particles. Clearly, the "colloidal" and the "atomic" results for D* collapse onto the same 
curve, which we denote by D* HS {(j)). One immediate and important consequence of this 



comparison is, for example, that Lowen's dynamic criterion for freezing 35( now applies 



for both, the atomic and the colloidal hard sphere liquid, i.e., the condition D* HS ((f)) ps 0.1 
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FIG. 3: Long-time self-diffusion coefficient Dl{(J)) of the hard-sphere fluid determined by molecular 



dynamics simulations 
normalized as D*((b) E 
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32l | . expressed in "atomic units" a^sT/M) 1 ' 2 (empty diamonds), and 



D L ((I))/D (0), with D°(» given by Eq. (I3~T3|) (full diamonds). The other 



nil symbols are the Brownian dynamics simulation results for D* from Refs. 
34J (circles). 



33] (triangles) and 



occurs at 



df) 
'hs 



0.494 in both cases. The comparison in this figure, however, is only 



one particular manifestation of the more general long-time dynamic scaling suggested by the 
present work, whose applications were also illustrated by the other results presented in this 
section. 

V. SUMMARY AND DISCUSSION. 

In this paper we have discussed the relationship between the dynamics of atomic and 
Brownian liquids, by describing the macroscopic dynamics of both kinds of systems within 
the same theoretical formalism. We have based this discussion on the application of the 
generalized Langevin equation formalism to the derivation of general memory-function ex- 
pressions for the (collective and self) intermediate scattering functions of an atomic liquid. 
The actual derivation, however, consisted in the review of the previous derivation 14J of 
the time-evolution equations for F(k,t) and Fs(k,t) of the corresponding Brownian fluid, 
keeping track of the effects of the solvent friction. At the end of such derivation the zero- 
friction limit was taken, to obtain the corresponding time-evolution equations for F(k, t) 
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and F s {k,t) of our atomic liquid (Eqs. (l2~29j) and (I2^0]l ). 

We then analyzed the long-time limit of these results for F(k, t) and Fg(k, t). The compar- 
ison of such overdamped expressions with the corresponding results in the case of a colloidal 
liquid, revealed the remarkable formal identity between the long-time expressions for F(k, t) 
and Fg(k,t) of atomic and colloidal liquids. As discussed in Sect. Ill, the fundamental dif- 
ference between these two cases lies in the definition of the diffusion coefficient D°; in atomic 
liquids it depends on temperature and density, and is given by the kinetic-theoretical result 
in Eq. (13.131) . whereas in colloidal liquids it is a constant, given by the density- independent 
Einstein- Stokes value in the absence of hydrodynamic interactions. Let us mention that this 
dynamic equivalence can also be inferred by the derivation of the (generalized) Langevin 
equation that describes the motion of representative tagged particles in an atomic liquid 



361 ] . The atomic-to-Brownian long-time dynamic equivalence thus seems to be a very ro- 
bust prediction, with important physical consequences. It implies, for example, that in an 
atomic system, the self-diffusion coefficient D° determined by kinetic theory plays the same 
role as the short-time self-diffusion coefficient D^ s ' in colloidal liquids. It also implies that 
the long-time dynamic properties of an atomic liquid will coincide with the corresponding 
properties of a colloidal system with the same S(k), provided that the time is scaled as D°t, 
with the respective meaning and definition of D°. 

In section IV we tested this observation by comparing the simulation results for Fg(k,t) 
obtained by both, molecular dynamics and Brownian dynamics, for the hard sphere system. 
As mentioned at the end of the previous section, one important consequence is that Lowen's 
dynamic criterion for freezing [35] now applies for both, the atomic and the colloidal hard 
sphere liquid. This result, taken together with the dynamic equivalence between soft- and 



hard-sphere liquids recently discussed in Ref. 26], further extends the application of this 
criterion to soft-sphere molecular liquids. The most relevant implications of this dynamic 
equivalence have been corroborated by the systematic comparisons between molecular and 
Brownian dynamics simulations of the sort illustrated in t his p aper. A summary of this 



analysis has been advanced in a recent brief communication [37 ]. 

We should mention, in addition, that in reality the validity of the present dynamic cor- 
respondence between atomic and colloidal liquids should extend over to colloidal systems 
involving hydrodynamic interactions, provided that the corresponding effects enter only 
through the value of the short-time self- diffusion coefficient D^(<p), which should then 
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play the role of a density- dependent D°, as suggested in [38j. Besides analyzing further 
these important predictions, we are in the process of applying the general expressions for 
F(k,t) and Fs(k,t) for an atomic liquid derived in this paper, to the development of a 
self-consistent scheme to calculate these properties. The intention is to extend to atomic 
he self-consistent generalized Langevin equation (SCGLE) theory of colloid dynam- 
161 ]. including the description of dynamic arrest phenomena |17H19|] and the recently 



liquids 

ics 15. 



developed first-principles theory of equilibration and aging 20|, |2l|. This, however, will be 
reported separately. 
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